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Abstract. We present the results of a narrow-band infrared imaging survey of a narrow strip (12' wide) around the 
galactic equator between 6° and 21° of galactic longitude aimed at detecting field stars with strong CO absorption, 
mainly late-type giants and supergiants. Our observations include follow-up low resolution spectroscopy (R = 980) 
of 191 selected candidates in the H and K bands. Most of these objects have photometric and spectroscopic 
characteristics consistent with them being red giants, and some display broad, strong absorption wings due to 
water vapor absorption between the H and K bands. We also identify in our sample 18 good supergiant candidates 
characterized by their lack of noticeable water absorption, strong CO bands in the H and K windows, and HKs 
photometry suggestive of high intrinsic luminosity and extinctions reaching up to Ay — 40 mag. Another 9 
additional candidates share the same features except for weak H2O absorption, which is also observed among 
some M supergiants in the solar neighbourhood. Interesting differences are noticed when comparing our stars 
to a local sample of late-type giants and supergiants, as well as to a sample of red giants in globular clusters 
of moderately subsolar metallicity and to a sample of bulge stars. A large fraction of the stars in our sample 
have Cal and Nal features markedly stronger than those typical in the local reference sample (both giants and 
supergiants), whereas the equivalent widths of the CO bands are similar or weaker. In this regard, our stars in 
the inner Milky Way disk display differences very similar to those identified by other authors between cool giants 
and supergiants near the galactic center and their counterparts in the solar neighbourhood. We propose that 
the systematic spectroscopic differences of our inner Galaxy stars are due to their higher metallicities that cause 
deeper mixing in their mantles, resulting in lower surface abundances of C and O and higher abundances of CN, 
which contribute to the strength of the Cal and Nal features at low resolution. Our results stress the limitations 
of using local stars as templates for the study of composite cool stellar populations such as central starbursts in 
galaxies. 

Key words. Stars: late-type, supergiants: Galaxy: abundances, disk, stellar contents, structure. 



1. Introduction 

As compared to the large amount of information avail- 
able nowadays on the inner regions of nearby galaxies, 
the Population I component of the central kiloparsecs 
of the Milky Way is still relatively unknown. Visual ex- 
Send offprint requests to: F. Comeron 

* Based on observations collected at the Centro Astronomico 
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Southern Observatory, La Silla, Chile (programme 71. B- 
0274(A)) 



tinctions towards objects located in the galactic plane 
only a few kiloparsecs from the galactic center, typically 
reaching several tens of magnitudes, restrict the obser- 
vations to the infrared and radio domains, or to X or 
gamma rays for objects associated with high energy phe- 
nomena. Moreover, a variety of factors such as distance, 
extinction, field crowding, line-of-sight confusion, or ab- 
sence of distinctive properties in the accessible wavelength 
range contribute to make difficult the identification or the 
detailed observation of individual objects. Observational 
programs such as large area surveys probing radio con- 
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tinuum and recombination line emission have supplied 
much of our present knowledge on the gas photoionized 
by massive stars (Lockman 1989, Hclfand et al. 1992 
Whiteoak lTMl Becker et al. HMfl Kuchar k Bania lTM21 
McClure-Griffiths ISUUTl Kolpak et al. |2TM|> . and on the 
nonthermal emission due to recent supernovae in the inner 
Galaxy (Leahy & Wu lB^l Green IllMllMill . Molecular 
line emission surveys have provided an overall picture of 
the distribution and physical conditions of giant molec- 
ular clouds (Dame et al. AMR IMiTl Chiar et al. IT5M1 
Bronfman et al. 11391)1 I2TJUU1 McQuinn et al. E^UUSjl . while 
mid-infrared surveys such as that carried out by IRAS 
have revealed a wealth of star forming sites (Wood k 
Churchwell 113591 Hughes k MacLeod IT3591 Codella et 
al.HUHl Comeron k TorraHHSl Egan et al. HMgjl . Recent 
near-infrared surveys like DENIS (Ruphy et al. 119971 
Lopez-Corredoira et al. 2001, Van Loon et al. |2003), 
2 MASS (Alard l2T)TTTl Ojha I21FTT1 Lopez-C orredoira et 
al. fM?^ or TMGS (Lopez-Corredoira et al. ITMfl [217)11 
Picaud et al. 2002J) are yielding extensive lists of stellar 
objects suitable for studies of galactic structure. Valuable 
as all those observational approaches are, each one has 
its own limitations and needs to be complemented by the 
others in order to provide a global picture of the young 
component of the inner parts of our Galaxy. On the other 
hand, this young component is an excellent tracer of pro- 
cesses of primary importance concerning the structure and 
the composition of our Galaxy, such as the dynamical re- 
sponse of the disk to the central bar, the location of the 
spiral arms, the current star forming activity, or the radial 
dependency of the chemical composition, thus justifying 
the efforts leading to the characterization of its properties. 

In this paper we present a study on a particularly 
promising, yet relatively unexploited, Population I tracer 
in the inner Galaxy. Late-type supergiant stars have in- 
teresting properties that make them particularly suitable 
for the study of aspects of the extreme Population I com- 
ponent that are either inaccessible or very elusive in other 
objects of similar age, such as the structure, kinemat- 
ics, and chemistry of the young galactic disk, while their 
numerical abundance can provide a measurement of its 
present-day massive star formation rate. They have ex- 
tremely bright magnitudes at infrared wavelengths, where 
the dust in the galactic disk is far less opaque than in 
the visible. As it turns out, a nearly complete census of 
late-type supergiants in our Galaxy is in principle a fea- 
sible task for a moderately large telescope equipped with 
a near-infrared camera, even if some of the targets may 
be as far as 25 kpc away and obscured by as much as 
50 mag at V. While the bright magnitudes of supergiants 
facilitate detection, the infrared spectral properties of cool 
photospheres facilitate their identification in conveniently 
designed infrared surveys. The ZVband spectra of objects 
with temperatures below ~ 3500 K are characterized by 
prominent CO bands appearing longwards of 2.29 /zm, 
thus allowing the selection of candidate late-type stars 
(red giant branch (RGB) stars, long period variables in the 
assymptotic giant branch (AGB), and red supergiants) by 



comparing narrow band images centered respectively on 
and off the CO band region. At low- and medium- reso- 
lution spectroscopy, many narrow atomic features in the 
H and K bands allow further investigation of their intrin- 
sic properties via their different sensitivities to tempera- 
ture, surface gravity, metallicity, and surface abundances. 
Moreover, infrared observations of cool, luminous stars in 
the inner galactic disk provide actual, individual stellar 
templates suitable for the interpretation of the integrated 
spectrum of systems undergoing massive star formation, 
such as central starbursts and the environments of AGNs 
(e.g. Origlia et al. HIM! Kotilainen et al. 11391)1 Origlia k 
Oliva 120001 Forster-Schreiber et al. 120011 Alonso-Herrero 
et al.HFjU, and are an important component in the mod- 
eling of the integrated light of entire systems (e.g. Fioc k 
Rocca-VolmerangeCnnZI Bruzual k Chariot EDOH- 

This paper describes a search for late-type luminous, 
obscured cool field stars in a 15°-long strip of the galac- 
tic plane of approximately 12' width in galactic latitude, 
consisting of narrow-band imaging and low resolution 
follow-up spectroscopy of the best candidates identified. 
Although our focus is on the identification criteria of cool 
supergiants we also have obtained as a byproduct a large 
number of spectra of cool RGB and AGB stars in the same 
general direction, which can be a useful resource for stud- 
ies of metallicity patterns across the inner galactic disk 
and for the construction of templates in integrated stellar 
populations spectra. 

The observations and data reduction process are de- 
scribed in Sectional The results are presented in Sectional 
where we discuss the properties of our sample and its com- 
parison to those of reference samples containing objects in 
the solar neighbourhood, in the galactic bulge, and in rela- 
tively metal-rich globular clusters. We discuss criteria that 
allow the definition of a sample of red supergiant candi- 
dates and present a list of the best ones identified among 
our observed objects. Section 01 discusses one of the main 
results of our work, concerning systematic spectral differ- 
ences that we find between our objects and those in the 
two reference samples. Section [S] summarizes our results. 

2. Observations 

2.1. Near-infrared imaging 

2.1.1. Observing strategy and filter choice 

The basis for our preliminary identification of possible 
cool supergiants consists of multiband images of narrow 
strips parallel to the galactic plane obtained with MAGIC, 
the near-infrared camera and spectrograph at the 1.23 m 
telescope at the German-Spanish Astronomical Center in 
Calar Alto. In imaging mode, MAGIC yields a pixel scale 
of 1" 15 per pixel resulting in a field of view of 4'9 x 4'9 
over the 256 x 256 pixel 2 array. Our imaging observations 
were obtained on two separate observing runs, carried out 
respectively on 3-13 July 1998 and 19-21 June 2001. 

The observations consisted of obtaining sequences of 
short exposures at positions separated by 1' intervals along 
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the galactic plane. Since the edges of the detector were 
aligned along the North-South and East- West directions 
and the series of images were obtained by moving the tele- 
scope along the galactic equator, stars lying at galactic lat- 
itudes near the edge of the coverage of the sequences were 
observed in fewer telescope pointings. At the beginning of 
each sequence the filter and the detector integration pa- 
rameters were set, and then the sequence was started to 
obtain a total of 180 overlapping fields covering an arc 
of 3° in galactic longitude. The starting positions corre- 
sponded to galactic longitudes of / = 6°, 9°, 12°, 15°, and 
18°, and galactic latitudes of b = 0', +4', and —4', thus al- 
lowing in principle overlap between contiguous sequences 
either in galactic longitude or latitude. Unfortunately, a 
major mechanical intervention on the mount of the 1.23m 
telescope between our two runs produced a slight misalign- 
ment of its polar axis and an inaccurate pointing model 
in our run of 2001. As a result, the sequences starting at 
b = —4' and the one at I = 15°, b = +4' arc not exactly 
parallel to the galactic equator and do not overlap over 
their entire lengths with contiguous sequences starting at 
the same galactic longitude. 

Four filters were used: a CO filter centered at 2.295 /im 
with a relative width AA/A = 1%, a i^-continuum narrow- 
band filter centered at 2.260 fim with AA/A = 2.6%, and 
the broad-band filters J (1.25 ^m, AA/A = 24%) and H 
(1.65 fim, AA/A = 18%). The choice of narrow-band fil- 
ters aims at allowing the detection of stars with large flux 
drops due to CO absorption, which is a characteristic fea- 
ture of late-type photospheres, by comparing the images 
taken through the narrow-band CO and X-continuum fil- 
ters (see below), while the broad-band niters were ini- 
tially intended to locate objects in color-color and color- 
magnitude diagrams. Exposure times per telescope posi- 
tion were 9 seconds in the CO filter; 6 in the i^-continuum 
filter; 4 in the H filter; and 8 in the J filter. These expo- 
sure times were obtained by stacking together as many 
individual exposures of 1 second each as needed. The ex- 
posure time per sky position is variable, as it depends on 
the number of telescope pointings containing that posi- 
tion as described above. It ranges from one to six times 
the exposure time per field, the latter value correspond- 
ing to stars lying close to the central galactic latitude of 
the sequence. Each frame was combined with the overlap- 
ping portion of the images taken before and after in the 
sequence in order to filter out remaining unmasked bad 
pixels or cosmic ray hits. 

In general preference was given to completing one se- 
quence in all niters before proceeding to the next one, in 
order to have observations as closely spaced as possible. 
However this had to be often overridden due to visibility 
constraints as the regions imaged, lying at declinations 
between —24° and —10°, are visible few hours per night 
from a mid-Northern site like Calar Alto. A consequence 
of this is that observations through the different filters are 
not consecutive in some cases, but separated by an interval 
of two days at most. Such time differences are nevertheless 
negligible by comparison to typical variability periods of 



evolved red stars, and we thus considered them entirely 
acceptable for the goals of our project. 

2.1.2. Data reduction and calibration 

The reduction of our imaging data, as well as of our spec- 
tra as described in Section 12.21 was carried out using 
IRAF 1 tasks and dedicated scripts. The images were re- 
duced using dome flat fields in each filter that were used 
also to construct bad pixel masks. Dark-subtracted, flat- 
fielded and masked images were sky-subtracted by con- 
structing a sky frame made out of the 20 frames taken 
closest in time, conveniently median-averaged and with 
the stellar images removed by sigma clipping. Next, the 
images obtained within each of the 6 = 0° sequences in 
the CO on-band filter, which we took as a reference, were 
mosaicked into a single large image. To this end, accurate 
offsets between consecutive pointings in a sequence were 
determined by automatically finding stars in the common 
area of overlapping images and using them as a reference. 
In case that a contiguous sequence (i.e., one starting at the 
same galactic longitude but a different latitude) was ob- 
tained during the same run, a similar procedure was used 
to align its component images with the reference 6 = 0° 
sequence: the position of each individual image along the 
parallel sequence with respect to the reference one was de- 
termined by using stars in the overlapping areas. This was 
not possible when sequences starting at the same galac- 
tic longitude were obtained on different runs, due to the 
pointing model inaccuracy outlined above. In such cases 
the sequence obtained on the second run was combined 
separately. The combined image in the CO filter was used 
as the reference for the alignment of each frame in the 
other filters, thus avoiding large misalignments between 
images in different filters due to accumulated small errors 
in the registering of consecutive images in each of them. 

Due to the crowdedness of the fields in the galactic 
plane and the undersampling of the point-spread func- 
tion, digital photometry with an undersized aperture of 2 
pixel (2" 3) radius was found to provide the best results. 
To compensate for the variability of the sky transparency 
conditions and for the variable image quality during the 
execution of any given sequence, instrumental zeropoint 
variations along the sequence were determined by com- 
paring the digital photometry on the brightest nonsatu- 
rated stars in the area common to each pair of overlap- 
ping frames. Zeropoints for each sequence in the J, H, 
and i^-continuum bands were estimated by comparing our 
instrumental photometry of bright nonsaturated stars to 
the 2MASS catalog, assuming that the narrow-band K- 
continuum magnitude can be approximated by the broad- 
band 2MASS K s magnitude. 

The outlined procedure for photometric calibration is 
appropriate to yield moderate accuracies in the J, H, and 

1 IRAF is distributed by NOAO, which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
under contract to the National Science Foundation 
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i-T-continuum bands that are sufficient for our purposes: a 
comparison of the instrumental and 2MASS magnitudes 
show typical scatters of ~ O.f mag around the average 
zeropoints, and a drift of the zeropoint of the same or- 
der at most along a given sequence. However, scatters 
of that order are too large for an appropriate measure- 
ment of the flux drop due to CO absorption, that we use 
as the primary criterion to select candidate cool stars in 
our sample. We thus used instead a self-calibrated method 
based on the fact that, in any random stellar field, most 
of the stars are expected to have a negligible drop in 
the flux at the position of the CO bands. We calculated 
local zeropoints of the photometric CO index (defined 
as [CO] phot = m(2.295) - m(2.26), where m(2.295) and 
m(2.26) are the magnitudes in the CO and if-continuum 
filters respectively) at each longitude by assuming that 
its value averaged over a large number of randomly se- 
lected stars must be very close to zero. This average was 
calculated considering the stars located within |AZ| = 10' 
of each longitude point and on the same galactic longi- 
tude band, and with weights inversely proportional to the 
square of the probable error in magnitude of each con- 
tributing star. The calculation was done in two steps, the 
first one including all the stars in the average and the sec- 
ond removing stars with a highly deviant [CO] p hot- I n this 
way, cool stars that could skew the distribution of [CO] p hot 
values could be selected as those stars standing out due to 
a value of that index greater than a given threshold above 
the average value of the surrounding stars. 

A master list of detected stars was obtained by 
coadding the images taken in all the filters and running 
DAOFIND ( Stetson ITWjl on the combined image. With 
this master list of stellar positions, we performed digital 
photometry with DAOPHOT on the images in each filter 
to obtain the multiband photometry of each star detected 
in the combined image. Many stars were not detected in 
some bands, either because they fell just outside the edges 
of the image due to small differences between the start- 
ing position of the sequence in each filter, because they 
were too faint in that filter to be detected, or because of 
confusion with a nearby source in the case of the images 
obtained under poorer seeing conditions. Finally, approx- 
imate stellar positions in J2000 coordinates were deter- 
mined by using as a reference the positions of unobscured 
stars in our images that could be identified in Digitized 
Sky Survey, and whose equatorial coordinates were avail- 
able from the USNO astrometric catalog. 

We must remark at this point that the conditions un- 
der which our observations were made resulted in an un- 
even depth of the catalog and other factors affecting its 
completeness. Seeing conditions during the execution of 
the observations ranged from ~ 1" to over 3", thus re- 
sulting in a regional dependence of the source confusion 
limits. Likewise, part of our sequences was affected by the 
presence of high altitude dust combined with the need for 
observing often at high airmasses, producing regional de- 
pendencies of the depth of the frames. The net effects of 
all these circumstances are rather uncertain and difficult 
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Fig. 1. Distribution of the magnitudes of the stars in our 
sample in the J, H, and if -continuum filters. Only stars 
with a < 0.3 mag in the corresponding band are consid- 
ered. The sharp cutoff at bright magnitudes indicates the 
saturation limit. 



to quantify, a fundamental limitation of the present survey 
that must be kept in mind. 
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Fig. 2. The distribution of the [C0] p hot index within three 
(H — K) color bins, showing that cool stars are abun- 
dant among the most reddened ones. Only stars with 
er([CO] p hot) < 0.1 have been plotted here. Note the use 
of a logarithmic scale in the vertical axis to enhance the 
wings of the distribution of indices, which is strongly 
peaked at [CO] p hot = 0. 

2.1.3. Properties of the photometrically selected 
sample 

Our final catalog of detected sources contains 435,445 en- 
tries corresponding to stars with their magnitudes de- 
termined to an accuracy better than 0.3 mag (excluding 
the contribution of the error in the zeropoint determina- 



tion) in at least one of the J, H and X-continuum filters. 
This level is typically reached at J = 15.5, H = 15.0, 
K = 13.0. The shallower limit at K is due to the use of 
a narrow-band filter to estimate the K magnitude. On 
the other hand, deeper limits in J and H are needed 
to estimate the reddening towards the heavily obscured 
stars that are the object of our program. Figure Ogives 
the histogram of source counts as a function of magni- 
tude, excluding stars with a > 0.3 mag in each band. The 
[COJphot can be usually measured down to an accuracy 
better than 0.1 mag for stars down to K ~ 10.0, as esti- 
mated from the scatter of the measured [CO] p hot values 
as a function of K . Although this is a rather modest limit 
in terms of depth, as we will see it is sufficient for our pri- 
mary purpose of detecting late-type, distant supergiants. 
Our catalog contains 23,388 stars with their [CO] p hot in- 
dices determined to that level of accuracy. These stars are 
plotted in Figure to illustrate the changing behavior of 
[COJphot with the (H — K) color index. The top panel 
shows a sample with relatively blue (H — K), which can 
be expected to be dominated by nearby stars with low 
reddening. The distribution of [CO] p hot is approximately 
symmetric around a central value of with a very minor 
fraction of the stars having an index significantly differ- 
ent from that value, in most cases due to measurement 
errors caused by field crowding or the presence of nearby 
companions. Such measurement errors can be reasonably 
expected to affect stars regardless of their (H — K) . The 
second panel, with moderately reddened stars, has a very 
similar appearance but shows hints of skewness towards 
high values of [CO] p hot- The skewness is more clear when 
we consider the reddest stars in our sample, where we ex- 
pect an increased contribution of luminous late-type stars 
that can be detected even at large distances and hence 
large amounts of foreground reddening. 



2.1.4. Sample selection for spectroscopic follow-up 

The imaging results described in the previous Section were 
then used to select the best late-type cool star candidates 
to be spectroscopically observed. This sample was com- 
posed of stars having [CO] p hot > 0.15 and a if-continuum 
magnitude brighter than 10. Our selection was further re- 
fined by visually inspecting the images in which the stars 
selected according to these criteria appear, and discarding 
stars whose photometry was uncertain due to proximity to 
the edge of the image (where contamination by unfiltered 
bad pixels or cosmic rays was more frequent) or confused 
by the proximity to nearby bright sources. Such selection 
produced a list of 204 targets to be spectroscopically ob- 
served, out of which three were discarded because of the 
presence of a close companion that was not resolved in 
the images from Calar Alto. The star numbers given in 
Tables ^ HI an d 0] refer to the entry number in this list, 
which is sorted by increasing right ascension. The posi- 
tions and broad-band photometry of the selected sources 
were obtained from the 2MASS Point Source Catalog, as 
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its astrometric and photometric accuracy is better than 
that obtained for our sample with the procedures de- 
scribed above. 

2.2. Near-infrared spectroscopy 

The follow-up spectroscopy of the best candidate cool 
stars selected as described above was obtained using the 
SOFI near-infrared camera and spectrograph at the ESO 
New Technology Telescope on La Silla, Chile, on the nights 
of 11, 12, and 13 July 2003. The setup and exposure pa- 
rameters used were identical for all the stars. We used the 
low-resolution red grism, which covers most of the H and 
K bands from 1.53 /im to 2.52 /im at a resolution R = 980 
with the 0" 6- wide slit that we used. Four spectra were ob- 
tained for each star following an AB/BA pattern, with 
each spectrum being the on-detector stack of three expo- 
sures of 10 seconds each. The B5V star HD 169033, which 
lies near the arc of the galactic plane where our targets 
are located, was observed at periodic intervals of approx- 
imately one hour to sample the evolution of the telluric 
absorption in the proximities of our field. On the second 
night of our run we observed both HD 169033 and the so- 
lar analog HD 172411 (spectral type G3V) consecutively 
and near their culmination, in order to minimize effects in 
the telluric features affecting both. 

For the reduction of the spectra, we started by can- 
celing out the sky emission lines and the detector back- 
ground by subtracting from each spectrum taken on the A 
position another taken on the B position, and viceversa. 
To correct for pixel-to-pixel variations we used a flat field 
obtained with the telescope enclosure closed, subtracting 
exposures taken with a continuum-emitting lamp respec- 
tively on and off. 

The wavelength calibration was carried out in two 
steps. We first obtained an exposure of a Xenon arc lamp, 
which produces a line spectrum conveniently sampling the 
wavelength interval covered by our observations, immedi- 
ately after observing one of our targets and keeping the 
instrument in the same position and setup, thus eliminat- 
ing any possible flexure effects or inaccuracies in reposi- 
tioning the slit. The arc lamp spectrum was wavelength- 
calibrated and the dispersion solution was applied to the 
sky spectrum of the object. In this way precise wavelenghts 
of the sky lines and line blends at the resolution of our 
observations were determined. The master calibrated sky 
spectrum was then used to wavelength-calibrate the rest 
of the observations. Each individual object spectrum was 
extracted, and the same aperture trace was used to extract 
the corresponding sky spectrum from an opposite frame 
in the AB/BA cycle. The lines identified in the master sky 
spectrum were then reidentified in the extracted sky spec- 
trum at the position of the object and used to determine a 
new dispersion solution, which was then applied to the ob- 
ject spectrum. Finally, the four individual sky-subtracted, 
flat-ficldcd, wavelength-calibrated spectra of each object 
were combined into a single spectrum. An identical pro- 



cedure was used for the telluric calibration star and the 
solar analog observations. 

To remove telluric features, the spectrum of each pro- 
gram star was ratioed by the spectra of each of the telluric 
calibration stars observed on the same night. The residual 
left after ratioing at the position of the strong CO2 fea- 
tures in the 2.00-2.05 ^m interval was used as a quality 
criterion to decide on the most suitable telluric calibra- 
tion star to correct for the telluric features in each pro- 
gram star. Since telluric features generally depend on both 
time and airmass the selected telluric calibration observa- 
tion was often, but not always, the one obtained closest 
in time to the observation of the program star. 

The last step was relative flux calibration and removal 
of photosphcric features in the spectrum of the star used 
for telluric correction. This was carried out by ratioing the 
consecutive observations of the telluric calibrator and the 
solar analog described earlier in this Section, and multi- 
plying this ratio by the flux-calibrated solar spectrum. The 
solar spectrum of Livingston & Wallace (1 991[1 was used 
for this purpose, with its resolution degraded to match 
that of our observations, following the procedure described 
by Maiolino et al. (1992|). The product of the telluric-over- 
solar analog ratio and the solar spectrum thus produced 
a multiplicative factor as a function of wavelength to be 
applied to our telluric-calibrated program objects in order 
to obtain their flux-calibrated spectrum in arbitrary units. 

Of the 201 stars observed, 194 turned out to display ob- 
vious CO absorption bands characteristic of spectral types 
K or M. We discarded the 7 spectra of stars photometri- 
cally misclassified as cool stars, as well as 3 additional 
spectra of faint stars showing hints of CO absorption but 
having an insufficient signal-to-noise ratio, which either 
prevented the determination of the equivalent widths dis- 
cussed in this paper or produced measurement errors not 
comparable to those of the other stars in the sample. 
Tabled (available in full electronically only 2 ) contains the 
entire sample of 191 stars presented in this paper listing 
their positions, 2MASS photometry, and measured equiv- 
alent widths of selected features (see Section I3.1JI . The 
numbering of the objects discussed in this paper refers to 
the entry number in that Table. The complete collection 
of reduced spectra used in the forthcoming discussions is 
also available electronically. 

As an additional step, we computed radial velocities 
for all our program stars by cross-correlating their iC-band 
spectrum to a zero-velocity cool star template spectrum 
obtained by combining the spectra of M-type giant and 
supergiant stars in the atlas of Wallace & Hinkle (1997). 
Doppler-corrected spectra were produced in this way for 
all our program stars, as well as for the stars in the Langon 
& Woods (2000) atlas that we use extensively as a refer- 
ence in the present work (see Section l3~2"|) . The Doppler- 
corrected spectra were the ones used for the derivation of 
the equivalent widths discussed in this paper. The accu- 

2 These data in electronic form can be accessed from the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 
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Table 1. Position, 2MASS photometry, and equivalent widths of selected features of all the observed stars 



Equivalent widths 

Number RA dec K s J - H H - K s CO(6,3) Nal Cal CO(2,0) I(H 2 0) 

(2000) (2000) (2MASS) (A) (A) (A) (A) (mag) 



001 


17:58:58.83 


-23:43:03.7 


10.36 


2.11 


0.87 


5.65 


4.34 


3.71 


15.64 


-0.185 


002 


17:59:17.61 


-23:48:02.3 


9.93 


2.09 


0.92 


8.34 


4.04 


3.25 


21.00 


-0.131 


003 


17:59:19.99 


-23:49:34.8 


9.86 


1.23 


0.53 


2.34 


2.31 


2.70 


11.33 


-0.100 


004 


17:59:20.46 


-23:51:11.7 


9.53 


2.47 


1.14 


9.45 


5.11 


3.42 


21.99 


-0.167 


005 


17:59:20.63 


-23:34:01.1 


9.64 


2.21 


0.95 


7.84 


4.32 


3.59 


20.44 


-0.180 



Note: the full version of this Table is available from the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 



racy of the radial velocities obtained is limited given the 
resolution of our spectra, and we estimate it to be approx- 
imately 30-40 km s _1 from the residuals of the dispersion 
solutions and the width of the cross-correlation peaks. 

3. Results 

3.1. Equivalent widths and spectrophotometry indices 

At the relatively low resolution available, the main fea- 
tures that are visible in our if-band spectra are the CO 
bandheads redwards of 2.293 /im, as well as the Nal fea- 
ture at 2.208 /im and the Cal feature at 2.264 //m 3 . Other 
atomic lines such as those of Fel and Mgl are also eas- 
ily visible, as well as weaker features due to CN. More 
features appear in the H band. The reader is referred to 
Kleinmann & Hall fiWty . Wallace & Hinkle ifTM?) and 
Hinklc ct al. ( 1995) for detailed line identifications in these 
two bands. The overall shape of the continuum is strongly 
affected by the presence or absence of the broad wings of 
the water absorption band centered around 1.9 yum, where 
the separation between the H and K bands lies. Figure |3 
shows some selected spectra representative of the variety 
found in our sample. 

A considerable number of slightly different definitions 
of spectral passbands useful for the quantification of dif- 
ferent features can be found in the literature (e.g. Lancon 
& Rocca-Volmerange 1 1 0921 Ramirez et al. 119971 Lancon 
& Wood lSUUUl Forster-SchreiberEIini Frogel et al. EMTl) . 
which are normally defined in the most convenient way 
given the spectral resolution available. Our choice has 
been to adopt the if -band indices defined by Frogel et 
al. l(2TM|l for the Cal, Nal, and CO(2,0) features, which 
are also directly comparable to those used by Ramirez et 

3 These features receive their names from the atomic species 
that normally dominate the absorption at those wavelengths in 
spectra with resolutions R ~ 1000 or lower. However, higher 
resolution observations reveal contributions also by CN and 
atomic species that can become important below 4,000 K 
(Ramirez et al. I1997(l . as discussed in Section 2] Although in 
this paper we refer to these features by their usual denomina- 
tion in the literature (Nal, Cal), the reader should be aware 
that these species may not be the only ones contributing to 
their strength. 



al. <|1997ll20~00b|l . They thus allow a straightforward com- 
parison between the spectral characteristics of our objects 
and those of the extensive samples of RGB stars observed 
by those authors in the solar neighbourhood (Ramirez et 
al. 11997(1 . the galactic bulge (Ramirez et al. 12000b}) . and 
globular clusters (Frogel et al. I2001f) (see Section 13.2(1 . 
Moreover, we have defined additional passbands to mea- 
sure Sil and CO(6,3) features in the H band, and Fel, 
Mgl in the K band. In all cases except for the i-T-band 
CO features the equivalent widths are measured with re- 
spect to the local continuum level interpolated between 
two featureless points on both sides of the feature. The 
lack of a continuum longwards of 2.29 /jm in our spectra 
forces us to use an extrapolation instead using two contin- 
uum points at shorter wavelengths. The definition of the 
on- and off-feature bands used in this paper is given in 
Table |2 Note that the interpolation or extrapolation us- 
ing two continuum bands makes the measured equivalent 
widths virtually insensitive to reddening, which is not the 
case for the [CO] p hot defined earlier. Typical uncertainties 
in the equivalent widths that we measure in our spectra 
are ±0.5 A, except for the CO(2,0) band where we esti- 
mate ±1.0 A due to the need to extrapolate from the bluer 
continuum. 

The width of the wings of the H2O feature separating 
the H and K bands requires a definition different from 
the equivalent widths discussed above to quantify their 
strength. We have thus defined a reddening-free water in- 
dex, /(H2O), based on the comparison of the fluxes in two 
bands located on the blue and red wing of the feature and 
the approximate continuum level measured on unaffected 
parts of the H and K bands. The locations and widths 
of these bands are also listed in Tabled where fjj, 
denote the fluxes on the if-band (blue) and iT-band (red) 
wing of the feature, and /jj, are the respective contin- 
uum fluxes. Figure0]shows the position of the bands used 
to measure /(H2O). Assuming that the extinction has a 
wavelenght dependence of the form A\ oc A -17 (Mathis 
HUSH), the quantity 



/(H 2 0) = -2.5 log ig- 2.081 log ^§ (1) 
in Jk 
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Table 2. Feature passbands and their reference continua 3.2. The reference samples 



Feature 



line/band 
A AA 

(H H 



blue cont. 
A AA 
(pm) H 



red cont 
Ao 
(Mm) 



"Low- and medium-resolution atlases of cool stars in the 



Sil 


1.5898 


0.0035 


1.587 


0.002 


1.5953 


CO(6,3) 


1.620 


0.0117 


1.6140 


0.006 


1.6325 


CN 


2.1335 


0.005 


2.127 


0.008 


2.139 


Nal 


2.2075 


0.007 


2.194 


0.006 


2.215 


Fel 


2.2275 


0.006 


2.2145 


0.006 


2.2327 


Cal 


2.2635 


0.011 


2.2505 


0.011 


2.271 


Mgl 


2.2807 


0.0035 


2.277 


0.004 


2.2865 


CO(2,0) 


2.2955 


0.013 


2.250 


0.016 


2.2875 



infrared have been presented in a number of studies (e.g. 
(fim) Kleinmann & Hall E38.6, Lancpn & Rocca-Volmerange 
Dallier et al. [TUMI Wallace & Hinkle [TUMI 115371 
Meyer et al. 119981 Forster-Schreiber 120001 Ivanov et 
al. 120041 see this latter reference for a more exhaustive 



Definition of passbands for measurement of H2 O band wings 



Ao 
(Mm) 



AA 

(Mm) 



fh 
fk 
f K 



1.675 
1.755 
2.075 
2.255 



0.05 
0.05 
0.05 
0.05 



Note: passbands for the measurement of the Nal, Cal, and 
CO (2,0) features are defined as in Frogel et al. §001) and 
Schultheis et al. (1277031 . 



O 




CD 



1.8 2 2.2 

Wavelength (microns) 

Fig. 4. The positions of the spectrophotometric bands de- 
fined in Section 3.1, plotted on representative spectra from 
our sample to show their usefulness in measuring the flux 
depression caused by strong water band wings. 



0.0025 
0.009 

004 °^ re l evan t work), with the general twofold purpose 
0045 °^ providing a tool both for the classification of individ- 
002 ua i objects and for the construction of templates of in- 
0.008 tegrated spectral energy distributions for population syn- 
0.007 thesis analyses. Among these atlases, a particularly suit- 
able one for the investigation of the properties of our sam- 
ple is provided by the extensive collection of cool stellar 
spectra in the visible and the near infrared published by 
Lancon & Wood PUUUI hereafter LW2000). Their sample 
provides a wide variety of spectra corresponding to differ- 
ent classes of luminous cool stars in different environments 
such as the solar neighbourhood, the galactic bulge, and 
the Magellanic Clouds, and includes multiepoch observa- 
tions for a large number of variable stars. Its resolution, 
R = 1100, closely matches that of the spectra presented 
in this paper, thus facilitating a direct comparison (see 
Figure [SJ. Classes represented in LW2000 include super- 
giants, thermally pulsating AGB stars of high and low 
amplitudes, non-pulsating RGB stars, and carbon stars. 

We have used the spectra from LW2000 for which 
those authors present observations in the H and K-b&nd 
and belonging to the samples of supergiants, oxygen- 
rich variables, bulge giants, and non-pulsating RGB stars. 
Unfortunately, few spectra in the H and K bands are 
available for the latter class of objects in the LW2000 
sample. This shortcoming can be largely circumvented 
by adding the collection of equivalent widths of the Nal, 
Cal, and CO(2,0) features in the JsT-band spectrum of 
bright giants measured by Ramirez et al. 0997) on spec- 
tra having a resolution similar to ours. We have also 
used the lists of equivalent widths of the same features 
in Ramirez et al. f2 00Qb| for as sample of over 100 
bulge RGB stars, and of RGB stars in the most metal- 
rich globular clusters observed by Frogel et al. <|2001[) . 
namely NGC 5927 ([Fe/H] = -0.37), NGC 6440 ([Fe/H] = 
-0.34), NGC 6528 ([Fe/H] = -0.17), and NGC 6553 
([Fe/H] = —0.25). Low-metallicity objects such as those 
in the Magellanic Clouds sample of LW2000 are not ex- 
pected to be significantly represented among our inner 
galactic disk objects, and have not been considered here. 
Likewise, there are clearly no carbon-rich stars in our sam- 
ple (which would be easily distinguishable thanks to the 
strong C2 feature at 1.77 /im) and we have thus excluded 
the objects of that class in LW2000 atlas from the com- 
parison with our results. 



is thus a reddening-free measurement of the depression in 
flux due to water vapour absorption on both sides of the 
gap separating the H and K bands, which we express in 
magnitudes for convenience. 



3.3. The H2O vs. CO diagrams 

A first classification of a sample of cool luminous star can 
be established on the basis of the relative importance of 
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Fig. 5. A comparison between the H and K spectra of one 
supergiant from the LW2000 library of spectra (above) 
and one supergiant candidate from our sample (below), 
showing their similar resolution and quality. The compar- 
ison also shows the reality of most of the features that 
can be seen in our spectra. The difference in continuum 
slopes, especially in the H band, is due to the higher ex- 
tinction towards the star of our sample. The features listed 
in Tableware noted. 

the most prominent molecular features in the H and K 
bands due to water steam and carbon monoxide. CO is 
to a first approximation a temperature indicator with the 
depth of the bands increasing as temperature decreases. 
However, the depth of the CO bandhead at 2.293 /jm 
is also sensitive to the surface gravity and is often used 
as a discriminator between giants and supergiants at a 
given temperature (e.g. Kleinmann & Hall ll98o1 Forster- 
Schreiber l2"UUUl Ivanov et al. l2TjUl|) . The CO band inten- 
sity is also sensitive to the carbon abundance in the at- 
mospheres of oxygen-rich stars and is thus an indicator 
of the importance of the dredge-up processes in the post- 
main sequence evolution of the star (Carr et al. EPOO). 
On the other hand, strong H2O absorption is formed in 
the cool, dense outer layers of large amplitude Mira-type 
and semiregular variables possessing very extended at- 
mospheres with a complicated structure determined by 



0.8 
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o 
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0.2 
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10 



20 



30 



EW[CO(2,0)] (A) 



Fig. 6. £W[CO(2, 0)] vs. 7(H 2 0) diagram for all the stars 
in our sample (eight-pointed asterisks), compared to the 
different classes of cool stars of LW2000: non-pulsating 
red giants (triangles), supergiants (squares), bulge giants 
(crosses), and oxygen-rich variables, separated between 
those having strong (circles) and weak (three-pointed as- 
terisks) water absorption. 

shocks propagating through them (Bessell et al. 119891 and 
references therein), and can be used as a discriminator 
between pulsating AGB stars and supergiants (Langon & 
Rocca-Volmerange ll992[l . both of which have overlapping 
luminosity ranges. 

Figurc[5]shows the distributions of our sample and the 
LW2000 one in the CO(2,0)-H 2 O diagram. The apparent 
overabundance of objects with strong water absorption in 
LW2000 is largely due to the fact that we have plotted 
in the diagram every single observation of each object in 
LW2000 as a separate point; therefore, any given star ob- 
served at different phases by LW2000 appears multiple 
times in the diagram. Figure [5] shows that our sample is 
strongly dominated by objects with the spectral character- 
istics of non-pulsating stars that occupy the locus of red 
giants and supergiants. The small scatter of /(H2O) for 
objects spanning a wide range of extinctions confirms that 
the extinction law chosen to derive Eq. (1) yields indeed a 
reddening-free index, whereas its steady decline between 
EW [CO (2, 0)] = 10 A and 20 A is due to the progressive 
departure of the H and K band spectrum from a Rayleigh- 
Jeans shape as the temperature decreases and the depth 
of the CO bands increases. Most of our objects have mod- 
erate values of EW[CO(2,0)] that are generally well be- 
low those of the sample of LW2000 supergiants, a feature 
that we will consider in the coming Sections. Finally, the 
objects with strong water absorption have a preference 
towards high £W[CO(2,0)], unlike in the LW2000 sam- 
ple where large positive /(H2O) indices are present over 
the whole range of £W[CO(2,0)] considered. The trend 
in our sample is qualitatively similar to that observed in 
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Fig. 7. £W[C0(6, 3)] vs. I(H 2 0) diagram comparing the 
stars in our sample with those of LW2000. Symbols are as 
in Figure 

Figure 13 of Lancon & Rocca-Volmerange Cl9iJ2J) where 
their objects split along two branches (corresponding to 
RGB stars and supergiants) at the highest EW[CO(2, 0)], 
although we do not find the gap between both branches 
suggested by the results of those authors. Unfortunately, 
no information on the H2O depression is available in the 
data published by Ramirez et al. Ijl997|) given their wave- 
length coverage restricted to the red part of the K band. 

Figure is the analogous of Figure El now replacing 
the CO(2,0) band by the CO(6,3) band at 1.62 /zm. Again 
cool supergiants of the LW2000 sample occupy the region 
of highest CO absorption, now with a larger overlap with 
our sample. A difference with Figure is that objects 
with /(H2O) > now populate a broad but clearly de- 
fined band so that the objects with the stronger water ab- 
sorption tend to have the lowest £W[CO(6, 3)]. The same 
behaviour in this respect can be seen in the LW2000 sam- 
ple and ours. A closer examination shows that not only 
the CO(6,3) band, but all the features in the H band are 
weaker among the stars with strong water band wings, in- 
dicating the existence of significant veiling correlated with 
the strength of the water feature. Some degree of veiling 
in the K band as well may be responsible for the relatively 
weak Cal and Nal features of stars with /(H2O) > in 
Figures El and El 

3.4. The Nal and Cal features 

Diagrams combining absorption in atomic lines and molec- 
ular bands having different dependencies on temperature, 
surface gravity, and metallicity offer in principle impor- 
tant diagnostics on the intrinsic properties of samples of 
individual cool stars and of composite populations, and 
have been discussed in a number of studies (see Ivanov et 
al. 120041 for a recent, comprehensive review). A compari- 
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Fig. 8. EW[CaI] vs. EW[CO(2,0)} diagram comparing 
the stars in our sample with those of LW2000. Symbols 
are in general as in Figure El but now the stars of our 
sample with /(H2O) > are marked as filled circles. Note 
the offset between the mean distribution of stars of our 
sample and other reference samples. 



son between our sample in the inner Galaxy, the sample of 
LW2000, and the samples of RGB stars in the solar neigh- 
bourhood, the bulge, and globular clusters described in 
Section f3. 21 can provide insight on possible intrinsic differ- 
ences of both populations and the way they reflect them- 
selves on measurable spectral features. 

We have plotted in Figures |S1 and El the relationship 
between the Nal and Cal features and the CO(2,0) band- 
head. In particular, the combination of the strength of 
a feature essentially dependent on the temperature such 
as Cal and one that combines dependence on both tem- 
perature and surface gravity such as CO(2,0) has been 
proposed as a luminosity class discriminator (Forster- 
Schreiber[2000, Ivanov et al.E 



Figures |H1 and El do not show tight overall correlations 
for the LW2000 sample, mainly due to the abundance in 
it of AGB stars whose complex, time-dependent atmo- 
spheric structure is expected to blur them. However, local 
supergiants dominate in a well defined region of the dia- 
gram characterized by the strongest absorptions in both 
the Cal and Nal features and CO(2,0). 

The loose correlation between the Cal or Nal features 
and CO(2,0) is much more apparent when only RGB stars 
and supergiants are considered, as shown in Figures ITDl to 
1131 where our sample is compared to the local supergiants 
of LW2000 and the RGB stars in the solar neighbour- 
hood (Ramirez et al. I1997|l . the galactic bulge (Ramirez 
et al. 2000b), and the relatively metal-rich globular clus- 
ters of Frogel et al. (2003). Clear systematic difference 
are noted between the stars in our sample on one side, 
and the local RGB sample of Ramirez et al. OTT)97 ■ and 



F. Comeron et al.: A search for late-type supergiants in the inner regions of the Milky Way 



11 



35 



30 



O 

o 



OX ° O 
! X » X 



•o" 

x '-l ^ 



• 



20 -■■ 



1 5 



A O 



m a o***; * ****** ^* 



I 



o °^ * 



*A * 



o 



o 



EWfNal] (A) 

Fig. 9. £W[NaI] vs. EW[G0(2, 0)] diagram comparing 
the stars in our sample with those of LW2000. Symbols 
are as in Figure [SJ Like in Figure [HI there is a clear offset 
between the mean distribution of stars of our sample and 
that of other reference samples. Note that a large fraction 
of our stars reach £W[NaI] values higher than those of 
most local supergiants. 
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Fig. 10. Same as Figure|5J but now comparing our sample 
to that of RGB stars in the solar neighbourhood from 
Ramirez et al. 1119971 four-pointed stars) and in relatively 
metal-rich globular clusters from Frogel et al. ( 2001) (open 
circles). The LW2000 sample of supergiants (open squares 
as in previous figures) is also included. 
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Fig. 11. Same as Figure^ but now comparing our sample 
to that of RGB stars in the solar neighbourhood from 
Ramirez et al. 1)19971 four-pointed stars) and in relatively 
metal-rich globular clusters from Frogel et al. I|2001|l 



(open 

circles) . The LW2000 sample of supergiants (open squares 
as in previous figures) is also included. Note that most of 
our stars have £W[NaI] greater than those of the coolest 
RGB globular cluster stars. 
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Fig. 12. Same as Figure|SJ but now comparing our sample 
to that of RGB stars in the galactic bulge from Ramirez 
et al. (|2000bl four-pointed stars). The LW2000 sample of 
supergiants (open squares as in previous figures) is also 
included. 



the sample of RGB stars in globular clusters of Frogel 
et al. H2001f) on the other, in the sense that the equiva- 
lent widths of the atomic features in our sample extend 
to higher values. Local RGB stars occupy essentially the 
same region of the diagrams presented in Figures an( i 



1111 as RGB stars in globular clusters, which are also plot- 
ted in those figures. This is especially clear in the Nal vs. 
CO(2,0) diagram, where the overlap in £W[NaI] between 
the giants and supergiants of the reference samples is mi- 
nor, and where only one RGB star belonging to NGC 6528 
(the highest metallicity cluster in Frogel et al. fiOOlJ) sam- 
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Fig. 13. Same as Figure^ but now comparing our sample 
to that of RGB stars in the galactic bulge from Ramirez 
et al. I|2000bl four-pointed stars). The LW2000 sample of 
supergiants (open squares as in previous figures) is also in- 
cluded. Note that most of our stars have £W[NaI] greater 
than those of the coolest RGB globular cluster stars. 

pie) appears with EW[NaI] > 4.0 A. The strengths of the 
CO(2,0) bandhead in our stars are otherwise similar to, 
or only slightly exceeding, the ones of the coolest RGB 
stars in the local and globular cluster samples. On the 
other hand, the galactic bulge sample has overall char- 
acteristics much closer to those of our objects, as shown 
in Figures ^] and ^] The largest measured equivalent 
widths of the CO(2,0) feature are similar to those of the 
local RGB star sample, which may be explained by their 
proximity to saturation and subsequent insensitivity to 
metallicity differences. However, RGB stars in the bulge 
do reach the high Cal and Nal equivalent widths that we 
measure in our sample and their loci in those diagrams has 
a much broader, although not complete, overlap with the 
one of our sample. The overlap with this sample is nev- 
ertheless poorer at low equivalent widths, where the local 
RGB stars from Ramirez et al. l)1997fl rather than bulge 
ones tend to populate the region occupied by the stars of 
our sample. This is to be expected since stars with lower 
equivalent widths are hotter, and therefore intrinsically 
fainter and more nearby on the average as indicated by 
the lower reddening of their colors, and can thus be con- 
sidered as an extension of the local sample. 

Metallicity differences between RGB stars in the so- 
lar neighbourhood and in the bulge are the most obvi- 
ous candidate to produce the systematic shifts between 
samples in the equivalent width diagrams. The quantifi- 
cation of these differences from low resolution spectra is 
however difficult, as the use of relationships calibrated on 
globular clusters implies their extrapolation outside their 
strict range of applicability (temperatures, luminosities, 
metallicities, feature strengths, star formation and chemi- 



cal pollution environments...) and even become dependent 
on the functional form adopted for the relationships, as il- 
lustrated by Figure 6 of Ramirez et al. (|2f)0f)bll . Moreover, 
the low resolution spectrum may be also determined by 
other effects that mimic metallicity effects, as discussed in 
Section O 

3.5. Color-magnitude diagram and luminosities 

The obvious differences in the intrinsic properties of the lo- 
cal sample of giants and the ones in the inner galactic disk, 
represented by both our objects and the bulge sample of 
Ramirez et al. ( 2000b ), leads us to be cautious about possi- 
ble differences also existing between local supergiants and 
those in the inner Galaxy that may render inapplicable the 
spectroscopic criteria allowing a separation between giants 
and supergiants in the solar neighbourhood. Indeed, such 
criteria do not seem to apply at least to giants and super- 
giants at the galactic center, as concluded by Schultheis 
et al. (2003) from their near-infrared study of a sample 
of cool, luminous stars detected in the ISOGAL survey 
(Omont 1996) with spectra obtained with the same in- 
strumentation as ours. Schultheis et al.'s identification of 
four very likely cool supergiants rests on their derived lu- 
minosities, which they estimate based on the assumption 
that all the sources in their sample lie at the distance 
of the galactic center. This is a reasonable working hy- 
pothesis in their case, given the high density of sources 
detected within a radius of less than 1°5 from the galac- 
tic center in which foreground contamination is expected 
to be negligible. Unfortunately we cannot make a similar 
assumption on the distance, and are thus forced to exam- 
ine more indirect means of estimating the luminosity of 
our stars in order to assess the possible presence of su- 
pergiants among them, without relying on spectroscopic 
criteria derived from the solar neighbourhood. 

Figure IT71 suggests that it is possible to make a rough 
luminosity estimate based on the 2MASS photometry 
alone. The figure plots the position of all objects in our 
sample in a (H — K$,Ks) color-magnitude diagram, in- 
cluding AGB stars characterized by their prominent steam 
bands (see Section 13 . 3|) . which are marked as full cir- 
cles. These stars are expected to be long-period variables 
roughly following the period-luminosity relationship for 
Miras and semiregular variables, which has been recently 
reexamined by Knapp et al. 1(277 03) using Hipparcos paral- 
laxes. Their absolute magnitudes range from Mk — —6.4 
for the shortest periods to Mk — —8.2 at the tip of 
the AGB, and thus overlap with the absolute magnitude 
range of cool supergiants, from which they can be oth- 
erwise easily separated spectroscopically. Concerning the 
brightest AGB stars, their absolute Mk magnitudes are 
expected to depend on the metallicity, primarily due to 
the cooler temperatures reached at the tip of the RGB 
by the most metal-rich populations as compared to oth- 
ers of lower metallicity, which result on a greater i^-band 
bolometric correction (Houdashelt et al. I2000[) combined 
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H-K s 

Fig. 14. Color-magnitude diagram for the stars in our 
sample. Filled circles represent stars with /(H2O) > 
(large-amplitude variables), expected to be intrinsically 
very luminous. Magnitudes are from 2MASS. The straight 
line represents the displacement due to an extinction of 
Ay = 10, using the Rieke & Lebofsky (1985) extinction 
law. 

with an essentially metallicity-independent bolometric lu- 
minosity (Sweigart et al.|T990). A moderate extrapolation 
of the metallicity-Mft- relationship established by Ferraro 
et al. (|2)00) from a sample of globular clusters to the so- 
lar metallicities indicates an expected Mk = — 7 for the 
brightest RGB stars, thus overlapping with the absolute 
magnitudes of the AGB stars with the shortest periods. 

The positions of AGB stars, whose absolute magni- 
tudes are confined to a range of ~ 2 mag, lie in a rel- 
atively well defined band of the (H — Kg, Kg) diagram 
whose slope results from the combined effect of extinction 
and distance on the apparent Kg magnitude. The fact that 
the scatter within this band is of the order of the scatter 
in the absolute magnitudes of Miras and semiregular vari- 
ables suggests that the amount of reddening in the H — Kg 
color can be used as a rough distance indicator. The intrin- 
sic colors of Mira variables are somewhat uncertain: Glass 
E595) finds noticeable differences in intrinsic colors 



et al. 

when comparing Miras in the Large Magellanic Cloud, the 
solar neighbourhood, and the galactic bulge. They find 
intrinsic H — K colors of bulge Miras (in the SAAO sys- 
tem) that range between 0.32 and 0.61, which is redder 
than the reddest RGB stars (Houdashelt et al. 120000 and 
supergiants (Ducati et al. 2001). Therefore, assuming the 
reddening in H — Kg as a proxy for the distance, RGB and 
supergiants having the same position in the Kg, H — Kg 
diagram as a given AGB star should be more distant, and 
therefore more luminous. 

The lower edge of the locus of our large-amplitude 
variables can be approximated by the straight line Kg = 
5.87 + 2.20(i7 — Kg), which allows us to define a reference 
magnitude Kq as 




H-K s 

Fig. 15. Color-color diagram for the stars in our sam- 
ple, based on 2MASS magnitudes. Filled circles represent 
stars with /(H2O) > (large-amplitude variables). Five- 
pointed stars represent our best supergiant candidates (see 
SectionEH). Note that the J magnitudes of stars with very 
red (H — Kg) colors are often unreliable due to the likeli- 
hood of contamination by nearby stars that are brighter in 
J. The straight line represents the displacement due to an 
extinction of Ay = 10, using the Rieke & Lebofsky (1985) 
extinction law. 



K =K S - 5.87 - 2.20(# - Kg) 



(2) 



whose positive or negative value indicates that the star is 
either below or above the lower edge of the AGB band, 
respectively. Although it may be tempting to directly re- 
late Kq to Mk, it must be noted that the straight line 
is just a convenient approximation to describe the locus 
of AGB stars in our sample and does not have a sound 
physical foundation. Indeed, it is easy to show that the 
assumption that Kq is simply the absolute magnitude off- 
set by a constant requires the distance modulus to be 
DM = 12.93 + 0A2E{H - Kg) (where E(H - Kg) is the 
amount of reddening, the coefficient 0.42 preceding it is 
derived from the Rieke & Lebofsky (1985) extinction law, 
and the constant after the equal sign follows from the as- 
sumption that the lower edge of the AGB band is defined 
by stars having Mk = —6.4), which does not allow for 
distances smaller than 3.8 kpc. Still, the use of Eq. (2) al- 
lows the split of our sample between " bright" and " faint" 
stars taking AGB stars as an appropriate dividing line. 

The fact that the bulk of stars in our sample falls below 
the Kq = line in Figure 1141 confirms that our sample 
is dominated by giants. However, there are some bright 
stars with different amounts of reddening occupying the 
AGB locus but with no or weak steam absorption; note in 
particular the group of faint, very red stars at (H — Kg) > 
2.0. These are some of the most promising M supergiant 
candidates in our sample. 
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The vast majority of the stars in our sample cluster 
along the narrow strip of the (H — K s ), (J — H) color-color 
diagram that traces the reddening vector (see Figure IT5|) . 
Some stars deviate from this strip apparently occupying 
the region of bluer-than- normal (J — H) colors at a given 
(H—Ks) that would indicate the existence of near infrared 
excess. These are however stars with very red colors whose 
faint J magnitudes are near or below the confusion limit 
for 2MASS in the galactic plane, for which the 2MASS 
J-band measurement is likely to be significantly contami- 
nated by nearby stars. We thus consider the J magnitudes 
as generally unreliable for very red objects, thus preclud- 
ing a discussion on the infrared excesses of different classes 
of objects. 



3.6. Equivalent widths and luminosities 

The use of Kq as defined in Eq. (2) allows us to study 
possible correlations between this quantity, purely derived 
from broad-band photometry, and the equivalent widths 
measured on the spectra. The most obvious correlations 
among both types of quantities are shown in Figures 1161 
and El While both figures are qualitatively similar in 
showing the increase of average luminosity with increas- 
ing CO band depth, there are some interesting differences 
between both. The CO(2,0) band depth vs. Kq diagram 
shows an approximately constant slope, with considerable 
scatter, over the whole range of equivalent widths cov- 
ered by our spectra. Fig. ll7l shows however a break in the 
trend around EW[CO(6, 3)] ~ 7 A with the distribution 
of points flaring up to brighter values of Kq. The clear dis- 
placement of the average distribution of stars with strong 
steam absorption towards lower EW[CO(6, 3)] is a direct 
consequence of the fact, already mentioned earlier (see also 
Figure 01 that strong steam absorption is correlated with 
weak CO(6,3) absorption, unlike what is observed in the 
case of the CO(2,0) bandheads. 

Less pronounced, but still obvious correlations with 
Kq are found for the Nal and Cal bands. We have also 
looked for possible systematic effects of ratios of equiva- 
lent widths with Kq that may provide spectroscopic dis- 
criminators, in particular the £W(CO(2, 0))/(EW(CaX) + 
£W(NaI)) ratio that Ramirez et al. (1997) proposed to 
discern between dwarfs and giants. The usefulness of that 
ratio in discriminating also between giants and supergiants 
is suggested by the results of Forster Schreiber |20iTQB , and 
has been investigated by Schultheis et al. ( 2003 ) for their 
galactic center sample. Unfortunately we can only confirm 
the negative results of Schultheis et al. (2003J), at least at 
the resolution common to our spectra, in which no hints 
of a systematic trend with the luminosity indicator Kq is 
seen. 

Other absorption features appearing in our spectra 
such as the Mgl (2.281 /xm) or Fel (2.228 ^m) lines may 
also be useful indicators of temperature, metallicity, or 
surface gravity. In particular we have attempted to use 
EW[CO(2,0)}/EW[MgI] as a possible surface gravity in- 
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Fig. 16. £W[CO(2,0)] vs. the Kq quantity defined in 
Eq. (2) for the stars in our sample. Filled circles indicate 
stars with I(H 2 0) > 0. 
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Fig. 17. Same as Figure but now with EW[CO(6,3)} 
in the horizontal axis. Stars with /(H2O) > tend to lie 
to the left due to the same veiling effect on the CO(6,3) 
band that causes the diagonal band in Figure 



dicator based on the results of Forster-Schreiber ((2000), 
which suggest two different characteristic values of that 
ratio for giants and supergiants at the coolest tempera- 
tures. Again we find no systematic trend in our results, 
although its possible existence may be hidden due to the 
weakness of the Mgl feature that becomes close to the 
typical ±0.5 A uncertainty in our equivalent width mea- 
surements, thus introducing a large random scatter in the 
ratios that we can determine. Similar uncertainties due to 
the measurement errors relative to the small equivalent 
width affect measurements involving the other lines, such 
as Fel (2.228 /xm) or Sil (1.590 /xm) 
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Fig. 18. A possible luminosity effect in the CN bands? 
In the top panel, spectra are sorted by the value of Kq 
(Eq. (2)), with the most negative values (implying on 
the average more luminous stars) at the top. The CN 
feature, marked by the dashed line, has a trend to be 
stronger in the most luminous spectra. The stars plot- 
ted here are the ones in our sample for which we measure 
24 < £W[CO(2, 0)] < 25, to remove possible temperature 
effects on the CN bands. However, such trend is not con- 
firmed for the general sample. The bottom panel shows 
the equivalent width vs. Kq, quantifying the effect hinted 
at by the top panel. 



LW2000 have noted that the strong CN features ap- 
pearing shortwards of 1 . 1 /im are characteristic of the su- 
pergiants of their sample. Unfortunately that spectral re- 
gion is too heavily obscured in the majority of our sources, 
thus preventing measurements of the strongest CN bands. 
However, there are multiple CN features appearing on the 
bluer side of the K band that may provide a useful al- 
ternative in heavily reddened objects. We find hints of a 
surface gravity dependence in the blend of CN lines near 
2.134 /im, as shown in Figure IT%1 Again, the weakness of 
the feature prevents its accurate measurement in many of 
our spectra, but higher resolution and S/N spectra might 
prove this to be an efficient way of producing clean sam- 
ples of giants or supergiants. We note that the compar- 
ison of CN and CO absorption bands provide measure- 
ments of surface abundances of CNO-processed elements 
bearing important implications on the internal structure 
and dredge-up history of the objects under study (Carr et 
al. 120001 see also Section 0} . 

3.7. Defining a sample of candidate supergiants 

As discussed in the previous Sections, we have failed to 
find in our spectra a discontinuity between the properties 
of giants and supergiants allowing a clear discrimination 
between both groups. Moreover, very few objects in our 
sample display simultaneously weak or absent H2O ab- 
sorption wings and very strong CO(2,0) absorption in the 
range covered by the supergiants of LW2000, thus sug- 
gesting that supergiants with spectroscopic characteris- 
tics closely matching those of the local sample are not 
represented in large numbers among the stars that we 
observed. Nevertheless, we do find many stars with Nal 
absorption similar or greater than those of local super- 
giants, no signs of H 2 absorption wings, moderate but 
not extreme CO(6,3) and CO(2,0) absorption, and photo- 
metric characteristics that suggest an intrinsic luminosity 
exceeding that of the brightest giants. These properties 
lead us to suspect that our sample does contain a number 
of cool supergiants in the inner Galaxy whose spectro- 
scopic characteristics present systematic differences with 
respect to local supergiants, being in fact closer to those of 
the galactic center supergiants, as noted in Sect ion l3~5l see 
also Section l4~T1 Taking together the indicators referred to 
above, we have produced a list of candidate supergiants 
in the inner Galaxy composed by the stars in our sample 
that fulfill the following criteria: 

— /(H2O) < —0.1, denoting the absence of signifi- 
cant H2O absorption wings. The negative threshold 
is quite restrictive, as weak steam absorption is com- 
mon among the coolest supergiants (Laneon & Rocca- 
Volmerange 19923 and the /(H 2 0) values of the typi- 
cal Mira-like spectra in our sample reach several tenths 
above zero. As a matter of fact, as many as seven of 
the LW2000 supergiant spectra would not be consid- 
ered as supergiant candidates on the grounds of this 
criterion (see Figure EJ. 
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— £W[NaI] > 4 A, thus selecting as supergiant candi- 
dates objects with a Nal feature stronger than that of 
the coolest giants in metal-rich globular clusters and 
in the solar neighbourhood. 

— £W[CO(2,0)] > 19 A and £W[CO(6,3)] > 7 A. We 
have already pointed out that these thresholds are well 
below the typical strengths of the CO features in the 
local supergiants, but they correspond to the values for 
which stars with Kq < appear in Figures fTTTl and IT7I 
As we will discuss in Section^] the possibility that the 
CO features are systematically weaker among bright 
cool stars in the inner Galaxy makes it advisable to 
use not-too-restrictive criteria on the CO features. 

— K < —0.1, to select stars whose broad-band photom- 
etry suggests high luminosity. The usefulness and lim- 
itations of K Q as a luminosity indicator have been dis- 
cussed in Section l3~51 and we stress again here that Kq 
is expected to be only loosely correlated with the ab- 
solute magnitude. We consider it however as a useful 
additional indicator of the possible supergiant nature 
of our candidates when used in combination with the 
other criteria described in this Section. The thresh- 
old chosen is intended to exclude stars that may have 
been placed above the Kq — line by inaccuracies in 
the photometry. 

The selection of the best cool supergiant candidates 
obtained through the application of these criteria thus 
contains 18 objects, listed in Table We have excluded 
one object, Star 191, due to its relatively blue (H — Kg) = 
0.40; as indicated in Section El the rough correlation be- 
tween broad-band colors and luminosity given by Eq. (2) 
must break down for blue colors indicative of low extinc- 
tion, and the characteristics of Star 191 are well compat- 
ible with those of a lightly reddened, nearby RGB star. 
It is interesting to note as a validation of these criteria 
that 7 out of the 8 reddest stars in our sample, which are 
heavily obscured stars but still bright enough to be in our 
sample and are thus good supergiant candidates based on 
their photometric properties alone, do fulfill all the spec- 
troscopic criteria as well. Likewise, Star Oil, which is the 
one that best matches the spectral characteristics of the 
LW2000 sample of M supergiants, is also included in our 
selection as a moderately reddened star and thus proba- 
bly not very distant, as may be expected from an object 
whose properties are similar to those of the local sample. 

As noted above, the criterion based on the almost com- 
plete absence of H2O absorption wings is deliberately con- 
servative, as the coolest local supergiants do show moder- 
ate absorption. We thus provide in Table 01 an additional 
list of M supergiant candidates, with a lower level of confi- 
dence than those provided in Table 01 with weak H2O ab- 
sorption. Such a table would contain most LW2000 local 
supergiants as well, although their EW[CO(2, 0)] would 
be greater on average. 

The Ks magnitudes of most objects in Tables GO and ^ 
and the reddenings implied by their H — K$ are consis- 
tent with distances to the Sun placing them well within 



the solar circle given the expected absolute K magni- 
tudes of cool supergiants. However this may not be true 
for the faintest, most reddened sources in those tables. 
For instance, assuming Mk = — 11 (which is very bright, 
although not unrealistic for a M supergiant; Elias et 
al.UHBSJ) and {H - K s ) = 0.3 for Star 072 would place it 
at 17.1 kpc from the galactic center on the opposite side of 
the solar circle, assuming a galactocentric distance of the 
Sun of 8 kpc, and obscured by Ay ~ 35 mag. Nevertheless, 
most of the other candidate supergiants are both brighter 
at K and less reddened, and the assumption that they 
are indeed located in the inner galactic disk can be ex- 
pected to be generally valid. This is even more so for the 
RGB stars of our sample, whose absolute magnitudes are 
fainter. 

4. Discussion 

4.1. Cool giants and supergiants in the galactic center 

Differences similar to those that we have discussed in the 
previous Sections between our sample of cool luminous 
stars and local samples have been already noted in pre- 
vious works on the galactic bulge, as already discussed in 
Section l3~4l as well as on the galactic center. The latter are 
of particular interest, since they include supergiants in ad- 
dition to RGB stars. Early infrared investigations of the 
stellar population at the galactic center revealed a high 
density of cool, luminous stars including supergiants and 
giants (Lebofsky et al. lT§g2l Sellgren et al. lWfjl . Similar 
to their hotter counterparts, the comparison of properties 
of cool massive stars in the galactic center with similar 
stars in the solar neighbourhood has been the subject of 
intense study to investigate how the widely different envi- 
ronment influences their formation and evolution. 

Sellgren et al. (I9i2| see also Blum et al. H996H already 
noted that a number of galactic center giants and super- 
giants display unusually strong Cal and Nal features when 
observed at low spectral resolution. Although differences 
between galactic center and solar neighbourhood luminous 
stars have been often interpreted as due to the effects of 
metallicity on stellar evolution, detailed abundance analy- 
ses have recently called this into question by showing that 
galactic center stars have metallicities very close to solar 
values (Carr et al. 120001 and references therein; Ramirez 
ct al. 12000a! Blum et al. l2llU5|) . 

Most interestingly in the context of our results is 
the detailed abundance analysis performed by Carr et 
al. (2000) on IRS 7, the brightest M supergiant in the 
galactic center. While spectroscopy at resolution similar 
to ours also shows strong Cal, Nal features for this star, its 
CO band at 2.293 /im is found to be significantly weaker 
than those of the local M supergiants a Ori and VV Cep, 
whose temperatures are very similar to that of IRS 7 as 
shown by model atmosphere fits. The study of Carr et 
al. H2000f) reveals a strong surface depletion of C and O and 
a corresponding enhancement of N with respect to local 
supergiants, while the combined abundance of C, N, and 
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Table 3. Candidate supergiants with no H2O absorption wings 
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Table 4. Candidate supergiants with weak H2O absorption wings 



Equivalent widths 

Number RA dec K s J — H H — Ks CO(6,3) Nal Cal CO(2,0) 7(H 2 0) 

(2000) (2000) (2MASS) (A) (A) (A) (A) (mag) 



064 


18:07:50.72 


-20:04:09.2 


7.77 


2.91 


1.58 


9.13 


4.21 


3.34 


19.58 


-0.063 


070 


18:08:37.35 


-19:50:05.3 


8.40 


3.24 


1.57 


9.05 


4.50 


3.68 


25.41 


-0.059 


074 


18:09:37.50 


-19:28:27.1 


10.07 


4.11 


2.03 


7.87 


5.16 


5.06 


26.27 


-0.098 


097 


18:14:53.10 


-17:00:51.3 


8.28 


2.75 


1.41 


11.40 


5.12 


3.90 


24.57 


-0.050 


121 


18:16:38.18 


-16:23:34.0 


9.09 


3.12 


1.62 


10.43 


5.76 


4.23 


25.73 


-0.024 


125 


18:16:55.97 


-16:09:54.3 


5.15 


1.91 


0.91 


9.46 


5.09 


4.59 


21.99 


-0.031 


152 


18:18:05.31 


-15:57:19.4 


6.48 


2.76 


1.44 


9.66 


5.59 


4.60 


24.29 


-0.035 


184 


18:25:11.20 


-12:28:40.2 


8.33 


3.02 


1.43 


10.04 


5.12 


4.40 


22.48 


-0.049 


194 


18:27:39.07 


-11:39:23.0 


7.09 


2.11 


1.03 


7.74 


4.59 


4.10 


23.88 


-0.045 



O adds to a value close to solar. Such altered abundances 
not only explain the weaker CO features (whose strength 
is essentially determined by the C abundance in O-rich at- 
mospheres) , but also the unusual strength of the Nal and 
Cal features in low resolution (R ~ 1000) spectra. Indeed, 
high resolution spectroscopy reveals that these atomic fea- 
tures are strongly contaminated by a large number of CN 
bands, which can even dominate the contribution to the 
feature in the coolest stars (Ramirez ct al. 1997). The en- 
hancement in the surface abundances of CNO-processed 
elements has the net effect of increasing the contribution 
of CN to the observed Nal and Cal features. 

The interpretation given by Carr et al. ( 2000 ) invokes 
deep mixing enhanced by fast stellar rotation as the ex- 
planation for the difference between galactic center and 



local supergiants. They also suggest that higher rotation 
velocities on the average, caused by the extreme prevailing 
conditions for star formation, could also explain the dif- 
ferences observed in other classes of high mass stars at the 
galactic center without requiring metallicities higher than 
solar. A similar explanation may also apply to the galactic 
center giants whose low-resolution spectral properties are 
also comparable. 

As noted in Scction fc.SI four more bona fide supergiant 
candidates have been recently identified among ISOGAL 
sources near the galactic center by Schulthcis et al. (J3003). 
Their classification as supergiants is based on the esti- 
mated absolute magnitude and the absence of water steam 
absorption. Three of these stars have EW[CO(2, 0)] in the 
20-24 A interval, well within the range covered by the su- 
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pergiant candidates in our sample and also by local RGB 
stars. The fourth one has £W[CO(2, 0)] < 28 A, which is 
also well below those of most local supergiants in LW2000 
sample. These results thus support the trend of galactic 
center supergiants to have CO bands weaker than those 
of local supergiants. 

4.2. Metallicity and systematic differences in feature 
strengths 

It is extremely unlikely that an explanation based on ro- 
tation such as the one proposed by Carr et al. (^000) can 
account also for the properties of our inner Galaxy sam- 
ple, since the formation conditions at distances of a few 
kiloparsecs from the galactic center should be much closer 
to those in the solar neighbourhood and no reason why 
our stars should be faster rotators is apparent. However, 
we note that the constraint of explaining the peculiar fea- 
tures of galactic center stars while not invoking supersolar 
metallicity most probably does not apply to our sample, 
which populates a region of the galactic disk where the ex- 
istence of a higher average metallicity is well established 
by a wide variety of tracers; see Chiappini et al. (2001 ) 
for an exhaustive summary of galactic abundance gradi- 
ent determinations for different elements. It is in fact very 
difficult to disentangle the effects of high metallicity from 
those of fast rotation on stellar evolution, and in particular 
on mixing in the interior of the star, as has been discussed 
by Schaerer lj3UUU|> . 

The resemblance between the distinctive properties of 
galactic center stars and those of our sample (greater 
strength of the Cal and Nal features, comparative weak- 
ness of the A"-band CO absorption relative to similar stars 
in the solar neighbourhood) leads us to consider it very 
plausible that the spectroscopic features of many of our 
stars can be explained also by the enhanced abundance 
of CNO-processed material in the stellar surface resulting 
from deep mixing, as suggested by Carr et al. ( 2000 ) for 
IRS 7 and possibly other galactic center objects. However, 
the fact that our stars most likely have a higher average 
metallicity than those of the solar neighbourhood and also 
those of the galactic center leads us to consider metallic- 
ity, rather than rotation, as the reason for the enhanced 
surface CNO element abundance. Higher metallicity im- 
plies higher opacity in the interior of the star, leading to 
a deeper penetration of the outer convective zone during 
the first dredge-up (Vandenberg & Smith IT$8"%1 Gratton et 
al.EOOO). In turn higher mass loss rates, enhanced by the 
higher metallicity, may also play a role by efficiently re- 
moving the outer envelope and rapidly make visible at the 
surface deep layers whose composition has been altered by 
CNO reactions (G. Meynet, priv. comm.) 

In view of our results and the above discussion, we 
refrain from deriving quantitative properties of our stars 
from the available spectroscopic material using extrapola- 
tions of relationships calibrated on populations consider- 
ably different from ours, such as the metallicity calibration 



of Frogel et al. (|20UT|) or the £W[CO(2, 0)]-ba sed temper- 
ature scale proposed by Ramirez et al. (1997). We find it 
necessary to defer a tentative determination of the intrin- 
sic properties of our sample (temperatures, metallicities, 
surface gravities, absolute luminosities...) to future work 
based on high resolution spectroscopy, which should allow 
quantitative analyses based on detailed atmosphere and 
stellar evolution models. 

5. Summary and conclusions 

In this paper we have presented the results of a search for 
cool, luminous stars in the inner region of the Milky Way 
(6° < I < 21°), selected by the signature of their strong 
CO bandheads longward of A = 2.293 /im in narrow- 
band imaging. Low resolution spectroscopy has been pre- 
sented for 191 stars confirmed to have moderate or strong 
CO bands and other signatures of cool spectra. Some of 
these stars are readily classified as large-amplitude vari- 
able candidates from the appearance of broad H2O ab- 
sorption wings. Our spectra also allow the measurement 
of a number of atomic and molecular features in the H 
and K bands. The addition of near- infrared photometry 
from 2MASS, which is available for all the stars in our 
sample, allows us to estimate rough luminosities. We find 
that our sample is dominated by red giant branch stars, 
but the photometric characteristics of some stars suggests 
luminosities in the range covered by the intrinsically more 
luminous large-amplitude variables, suggesting that they 
are red supergiants. The spectroscopic characteristics are 
consistent with such classification. 

We have compared the spectroscopic results obtained 
for our sample to those derived from the atlas of near- 
infrared spectra of cool luminous stars of Langon & 
Wood (2000), which is mostly composed of stars in the 
solar neighbourhood and bulge red giants, as well as to 
the samples of nearby red giants of Ramirez et al. f.1997) 
and bulge giants of Ramirez et al. ( 20003- We have also 
made a comparison to a sample of red giants in globular 
clusters with metallicities in the -0.37 < [Fe/H] < -0.17 
from Frogel et al. 1|2001J) . The results of this comparison 
may be summarized as follows: 

— The local sample of red supergiants displays spectra 
with CO bands significantly stronger than those of 
the stars in the inner Galaxy sample, especially the 
CO(2,0) band. The difference is less pronounced con- 
cerning the CO(6,3) band, where there is significant 
overlap between our sample and the local one. We find 
only one star in our sample whose £W[CO(2, 0)] is 
well in the range of those of local supergiants. 

— The Cal and Nal features in the K band are on 
the average as strong, or even stronger, than those 
of local supergiants. As compared to local and glob- 
ular cluster RGB stars, most stars in our sample 
have -EW^Cal] and, especially, £W[NaI] well exceeding 
those of the coolest stars in those samples. However, 
similarly strong Cal and Nal features are observed 
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among bulge giants. The net result of this effect and 
the one described in the previous item is a well defined 
shift of the band defined by our sample in the EW [Cal] 
vs. £W[CO(2,0)] or the £W[NaI] vs. £W[CO(2,0)j 
with respect to the average loci of the reference sam- 
ples of the solar neighbourhood. 

— An application of combined photometric and spec- 
troscopic criteria allows us to define a sample of 18 
good M supergiant candidates having no noticeable 
H2O absorption, strong CO absorption both in the 
H and K bands, and H, K$ photometry similar to 
that of large amplitude variables, suggesting that they 
are also intrinsically very luminous. Most of them 
are considerably reddened by extinctions reaching up 
Ay ~ 40 mag. Nine additional candidates with weak 
H 2 absorption wings are also identified. 

— A comparison between the spectroscopic characteris- 
tics at low resolution of our best supergiant candidates 
and of the rest of stars in our sample confirms results 
of other authors about the unsuitability of spectra at 
R ~ 1000 or lower to clearly discern between red giants 
and supergiants, complicated with the interpretation 
of the measurements of the CO and atomic features, 
and in particular the actual species that contribute to 
the latter. 

The discussion of our results (Section |3J is centered 
on the fact that the low-resolution characteristics noted 
in our sample and summarized above (strong Nal and 
Cal features, probably due to their contamination by CN 
lines, and comparatively weak CO absorption) are shared 
by giants and supergiants in the vicinity of the galactic 
center. We have referred to detailed abundance analyses 
by other authors of the galactic center stars, which in- 
dicate solar metallicities but altered surface abundances 
of CNO-processed elements with respect to local super- 
giants, resulting in depleted C and O and enhanced N. 
Systematic differences in rotational velocities, implying 
deeper mixing in galactic center stars, are invoked to ac- 
count for the different surface abundances while keeping 
a solar metallicity. While we find the ultimate reason of 
systematic differences in rotation velocities highly implau- 
sible to explain the characteristics of our sample, we still 
consider that deeper mixing and the resulting patterns in 
surface abundances are an attractive explanation. We pro- 
pose that higher metallicities in the inner regions of the 
Galaxy are responsible for the deep mixing, probably as a 
consequence of the increased interior opacity. 

Higher resolution spectroscopy in the infrared (the 
only spectral region accessible for the study of the photo- 
spheres of such heavily reddened objects) will be needed 
to perform detailed abundance analyses as well as to pos- 
sibly separate giants and supergiants in a more reliable 
way than was possible with the material currently avail- 
able. Nevertheless, interesting conclusions may already be 
drawn for the application of the results of our work. On 
the one hand, we confirm the feature, already suspected 
from studies of galactic center stars, that very strong CO 



absorption longwards of 2.293 /im is a distinctive property 
of local supergiants but may not apply to supergiants in 
other environments. On the other hand, the usefulness of 
local giants and supergiants as templates for population 
synthesis studies must be regarded with some caution tak- 
ing into account our results and the ones referred to on the 
galactic center. Indeed, common domains of application of 
population synthesis are central starbursts and the stellar 
environments of AGNs, where the conditions for star for- 
mation and evolution may be better represented by those 
in the inner regions or the center of the Milky Way. The 
identification of a large sample of cool, luminous field stars 
such as the one presented in this paper is a necessary first 
step for the understanding of their physical and chemical 
properties that may emerge from the combination of fu- 
ture high resolution observations and detailed modelling. 
Such work will provide a database useful to improve our 
knowledge on the star forming history and chemical evo- 
lution of the inner galactic disk in the proximity of the 
galactic plane, as well as of a variety of extragalactic envi- 
ronments where large scale massive star formation is tak- 
ing place. 
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Fig. 3. Representative examples of cool stellar spectra found in our sample. The interval between 1.80 /im and 2.05 /^m 
is contaminated by strong telluric absorption bands and has been clipped off. The spectra in the top left panel show 
relatively strong CO absorption bands, but no signs of depression due to the wings of the broad water steam feature 
centered near 1.9 /im. The top right panel shows a sample of stars with similarly strong CO, but now also with a 
conspicuous H2O feature that we classify as large amplitude variables following LW2000. The stars in the bottom set 
of spectra are drawn from the sample of candidate supergiant stars (see Section l3.7|) . The spectra are normalized to 
the flux between 2.20 and 2.25 /im, and the offsets between consecutive spectra are of 0.5 in these normalized flux 
units. The number accompanying each spectrum refers to the entry number in Tabled 



